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SECTION  I 


INTRODUCTION 

During  the  past  twenty  years  a  coordinated  program  of  research  on 
problems  related  to  the  fluid  dynamics  of  axial-flow  turbomachinery  has 
been  developed  at  Iowa  State  University  in  the  Turbomachinery  Components 
Research  Laboratory.  Numerous  projects  funded  by  the  National  Aeronautic 
and  Space  Administration  (NASA),  the  United  States  Air  Force  Office  of 
Scientific  Research  (USAF/AFOSR) ,  the  United  States  Air  Force  Aero 
Propulsion  Laboratory  (USAF/AFAPL) ,  the  Naval  Air  Systems  Command,  and 
other  government  and  industrial  organizations  have  been  completed. 
Beginning  1  October  1978,  the  Iowa  State  University  Department  of  Mechani 
cal  Engineering  faculty  involved  with  turbomachinery  research  were 
funded  through  AFOSR  Contract  F49620-79-C-0002  to  develop  a  multi-task 
program  concentrating  on  problems  in  the  aerodynamics  of  advanced  axial- 
flow  compressors,  fans,  and  turbines.  Five  projects  were  scheduled 
during  the  initial  two-year  contract  period. 

Task  I:  Analytical  and  Experimental  Investigation  of  Intra¬ 

passage  Flow  in  a  Constant  Mean  Radius  Rectangular 
Cross-Section  Passage  Representative  of  Passages  in 
Turbomachinery 

Task  II:  Evaluation  and  Management  of  Unsteady  Flow  Effects  in 

Axial-Flow  Compressors 

Task  III:  Development  of  a  Design-System  Oriented  Deviation 

Angle  Prediction  Equation  for  Advanced  Compressor  and 
Fan  Configurations 


Task  IV:  Experimental  Study  of  Merits  of  Specific  Flow  Path 

Geometry  Changes  in  Controlling  Secondary  Flows  in 
Axial-Flow  Compressors 

Task  V:  Definition  of  Experimental  Programs  and  Facilities 

Appropriate  for  University  Turbomachinery  Research 
Programs 

Section  II  of  this  report  reviews  the  progress  made  in  reaching 
the  objectives  in  each  Task. 
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SECTION  II 


RESULTS  OF  THE  RESEARCH  PROGRAM 

Each  of  the  five  projects  in  the  contract  program  is  discussed  in 
this  section.  Task  I  originated  in  internally-supported  work  at  Iowu 
State  University  and  is  based  on  earlier  analyses  carried  out  with  support 
from  the  Pratt  &  Whitney  Group  of  the  United  Technologies  Corporation. 

Task  II  is  a  continuation  of  research  funded  by  AFOSR  since  1975.  Task 
III  continues  computation  and  analysis  initiated  in  1974  with  funding 
from  NASA  and  extended  under  an  AFOSR  grant  in  1978.  Task  IV  evolved 
from  discussions  with  Dr.  Arthur  J.  Wennerstrom  of  AFAPL,  and  is  intended 
to  result  in  fundamental  design  data  that  will  be  usable  in  future  high- 
efficiency  multistage  axial-flow  compressor  development.  Task  V  is  an 
independent  study  of  the  current  capabilities  of  university  turbomachine 
research  groups  in  the  United  States,  and  of  project  areas  in  which 
these  groups  can  make  significant  technical  contributions. 


1.  TASK  I:  EXPERIMENTAL  AND  ANALYTICAL  INVESTIGATION  OF  INTRAPASSACE 
FLOW  IN  A  CONSTANT  MEAN  RADIUS  RECTANGULAR  CROSS-SECTION  PASSAGE 
REPRESENTATIVE  OF  PASSAGES  IN  TURBOMACHINERY 

1.1.  Introduction 

The  application  of  two-dimensional  or  blade-element  design  systems 
for  turbomachinery  cascades  is  often  subject  to  failure  in  cases  where 
large  secondary  flows  result  from  the  interaction  of  endwall  boundary 
layers  with  the  remainder  of  the  flow  field.  This  is  especially  true  in 
low  aspect  ratio  geometries  where  the  influence  of  the  endwalls  may  be 
felt  throughout  the  entire  flow  passage.  It  has  been  recognized  for 
some  time  that  extensive  additional  information  is  required  to  help  sup¬ 
plement  practical  cascade  design  systems,  with,  ultimately,  three- 
dimensional  flow  analyses  forming  the  basis  for  such  design  systems. 

For  this  reason,  a  detailed  investigation  of  intrapassage  flow  in  a 
large-scale,  curved  rectangular  cross  section  passage  representative  of 
turbomachinery  passages  was  undertaken.  The  experimental  results  thus 
obtained  may  provide  further  knowledge  of  the  complex  internal  flow 
structure  involved  and  may  also  be  used  in  evaluating  and  checking 
theoretical  analyses  as  they  are  developed  in  future  work. 

1.2.  Experimental  Facility  and  Passage  Test  Rig 

In  the  subject  curved  passage,  significant  secondary  flow  can  occur 
without  additional  complications  of  cascade  leading  and  trailing  edge 
effects  and  tip  clearance  flows.  Also  with  low-speed,  incompressible 
flow  testing  the  essential  details  of  the  three-dimensional  flow  field 
are  present.  In  the  present  investigation,  the  test  passage  consists  of 
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a  90-deg.  bend  of  constant  mean  radius  fo] lowed  bv  a  straight  discharge 
section  as  shown  in  the  schematic  diagram  in  Fig.  1.  The  actual 
test  passage,  inlet  and  discharge  sections,  and  associated  instrumenta¬ 
tion  are  shown  in  the  photographs  in  Figs.  2,  3  and  4. 

Air  flow  to  the  test  rig  is  supplied  by  the  Turbomachinery  Low-Speed 
Flow  Loop  which  consists  of  a  large  centrifugal  blower  with  design  flow 
rate  of  20,000  cfm  at  20  inches  of  water  head  rise,  variable-speed  drive 
and  motor,  a  piping  system  of  30-inch  diameter  pipe  including  metering 
section  and  supply  header,  and  a  large  diameter  plenum  discharge  section. 
Contained  in  the  plenum  is  a  series  of  screens  and  a  honeycomb  to  ensure 
satisfactory  flow  uniformity  and  turbulence  level.  The  flow  loop  is 
operated  open-circuit  for  the  present  testing  with  flow  exhausting  from 
the  test  rig  to  atmosphere. 

The  test  section  proper  consists  of  the  90-deg.  bend  and  discharge 
passage  as  formed  by  the  upper  and  lower  endwalls,  and  by  the  inner 
(suction)  and  outer  (pressure)  sidewalls.  As  can  be  seen  in  Fig.  2,  the 
upper  endwall,  which  was  fabricated  from  19-mm  thick  aluminum  tooling 
plate,  contains  12  traverse  slots  across  the  flow  passage  at  uniformly 
spaced  stations  in  the  flow  through  the  bend  and  straight  discharge 
sections.  Brass  tooling  bar  stock  (9.525  mm  by  19.050  mm)  was  used  for 
the  traverse  slide  bars  in  the  aluminum  endwall.  These  slide  bars, 
which  were  precision-fitted  in  the  endwall  to  produce  a  smooth  inner 
wall  for  the  passage,  were  used  to  carry  a  five-hole  directional  pressure 
probe  and  a  static  tap  fixture  for  mapping  out  the  details  of  the  flow 
in  each  traverse  plane.  The  slide  bars,  each  actually  composed  of  two 
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BOTTOM  PLEXIGLASS  ENDWALL 


side  of  bend.  Task  I. 


bars  riding  in  inner  and  outer  slots  in  the  endwall,  could  be  positioned 
across  the  passage  as  desired  and  clamped  to  seal  against  flow  leakage. 

The  bars  and  slots  necessarily  extended  beyond  the  sidewalls  of  the 
passage  in  order  to  keep  the  endwall  sealed  with  the  probe  positioned 
toward  either  side.  It  was  found  that,  although  the  slide  bars  fit 
tightly  in  the  endwall  slots  and  against  the  sidewalls,  they  could  be 
easily  positioned  manually  after  having  been  lubricated  with  a  lighL 
grease.  No  gasketing  of  the  slide  bars  against  the  slots  or  as  extra 
sealing  against  the  upper  edge  of  the  sidewalls  was  found  to  be  neces¬ 
sary.  The  lower  endwall  of  the  test  section  is  made  of  19-nun  thick 
plexiglass  plate  for  flow  visualization  and  is  attached  to  the  lower 
edge  of  the  sidewalls  by  quick-release  bolts  and  mounting  brackets  for 
easy  access  to  the  flow  passage.  Dimensional  details  of  the  test  section 
and  traverse  planes,  and  the  coordinate  scheme  used  for  flow  testing  art- 
given  in  Appendix  A. 

The  sidewalls  of  the  test  section  were  made  up  of  laminated  mahogany 
boards.  The  sidewalls,  which  can  be  seen  in  the  side  view  of  the  test 
section  in  Fig.  3,  or  in  the  view  of  the  test  section  discharge  in  Fig. 

4,  are  constant  height  of  40.6  cm  and  approximately  10  cm  thick.  Fabrica¬ 
tion  of  the  sidewalls  posed  a  difficult  machining  problem  because  of  the 
large  work  pieces  involved  and  limited  machine  and  milling  cutter  capa¬ 
city.  Machining  was  finally  carried  out  by  constructing  each  wall  from 
three  individually  machined  layer  segments.  In  forming  up  the  passage 
with  the  sidewall  inner  profiles  a  constant  mean  radius  was  maintained; 
however  to  help  accelerate  the  flow  through  the  bend,  a  slight  contraction 
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in  passage  width,  or  cross  sectional  area,  was  incorporated  through  the 
full  90  degrees  of  the  bend.  The  passage  width,  w,  is  given  by  the 
relation, 


—  -  1  -  —  (  L-a)  (4a-sm4't) 
w  T\ 

o 

where  w^  =  w(u  =  0),  and  a  is  the  (normalized)  assigned  width,  w(  <  -  •■/.<) 
With  this  relation,  the  transition  in  passage  width  between  the  bend  and 
the  constant  width  inlet  and  exit  sections  is  continuous  through  the 
second  derivative.  In  each  of  the  12  traverse  planes,  16  static  taps 
uniformly  spaced  in  each  sidewall  along  with  the  static  tap  fixture 
installed  in  the  slide  bar  were  used  to  obtain  static  pressure  distribu¬ 
tions.  The  back  side  of  these  taps  can  be  seen  in  the  side  view  of  the 
test  rig  in  Fig.  3.  The  static  taps  were  constructed  by  cementing  in 
place  short  segments  of  16  ga.  (1.2  mm  I.D.)  hypodermic  tubing  in  drilled 
holes  carefully  aligned  normal  to  the  sidewall  profile  at  the  traverse 
stations.  Care  was  taken  that  the  inserted  tubing  was  finished  flush  to 
the  sidewall  and  the  hole  deburred  to  obtain  a  clean,  sharp  static  tap. 

As  noted  in  the  schematic  of  the  test  rig  in  Fig.  1,  a  transition 
section  from  circular  to  rectangular  cross  section  followed  by  a  boundary 
layer  generator  section  and  adaptor  section  introduce  the  flow  into  the 
test  section.  Again,  these  sections,  which  were  constructed  from  molded 
fiberglass  and  plexiglass  plate,  can  be  seen  in  side  view  of  the  test  rig 
in  Fig.  3.  The  boundary  layer  generator  section  consists  of  a  removable 
screen  holder  segment  by  which  curved  and/or  slanted  uniform  mesh  screens 
can  be  installed  into  the  flow  without  disassembly  of  the  test  rig.  Thus 


controlled,  symmetric  two-dimensional  velocity  profiles  across  the  passage 
from  one  endwall  to  the  other  can  be  easily  produced  in  the  inlet  flow 


to  the  test  section.  The  following  adaptor  section  has  adjustable  side- 
wall  bleed  doors  to  provide  further  control  on  the  flow  into  the  test 
section.  Also  in  the  adaptor  section  a  pitot  probe  could  be  positioned 
to  the  center  of  the  cross  section  of  the  passage  to  determine  free- 
stream  dynamic  and  static  pressures,  q  and  pso,  respectively.  This  was 
done  when  setting  the  flow  rate;  the  pitot  was  of  course  removed  prior 
to  the  actual  testing.  Generally,  tests  were  run  at  a  plenum  gage 
pressure  of  approximately  235  mm  H0O,  or  corresponding  maximum  flow 
velocity  into  the  test  section  for  test  results  presented  later  of  about 
AO  m/s.  The  test  Reynolds  number  based  on  this  velocity  and  the  span  of 
the  passage  cross  section  is  1.02  x  10^. 

1.3.  Data  Aquisition  and  Experiment  Control  System 

As  a  major  effort  of  the  present  investigation,  an  automatic  data 
acquisition  and  experiment  control  svstem  was  set  up  for  tests  using  the 
large-scale,  low-speed  flow  loop  (see  Ref.  1  for  complete  description  of 
the  data  system).  This  system  was  designed  with  /ersatility  and  flexi¬ 
bility  in  mind  for  general  application  to  testirg  programs  in  the  labora¬ 
tory.  With  this  system,  experimental  data  are  recorded,  stored,  and 
instruments  controlled  automat ical Iv  as  programmed  in  the  test  procedure. 
In  addition,  data  are  collected  at  a  high  rate,  facilitating  on-line 
reduction  of  data  and  checking  of  experimental  results;  in  the  process, 
intervention  of  human  error  in  data  collection  and  reduction  is,  of 
course,  minimized. 
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A  schematic  layout  of  the  data  system  can  be  seen  in  Fig.  5.  Unknown 
pressures  supplied  from  the  pressure  probe,  pressure  reference  system,  or 
through  the  fluid  switch  as  indicated  are  fed  into  the  pressure  scanner 
where  a  pressure  transducer  converts  pressure  differences  into  voltage 
differences.  The  voltage  across  the  transducer  is  amplified,  filtered 
and  fed  into  the  voltage  scanner.  By  selecting  the  appropriate  scanner 
channel,  the  amplifier  output  is  connected  as  input  to  the  digital  volt¬ 
meter.  Selection  of  other  scanner  channels  provides  the  connection  of 
probe  position  voltages  to  the  voltmeter  or  the  activation  of  the  pres¬ 
sure  scanner  and  probe  positioner.  Once  the  amplifier  output  is  con¬ 
nected  to  the  voltmeter,  the  voltmeter  is  triggered  to  take  the  measure¬ 
ment.  The  reading  is  then  output  to  the  system  PET  controller/computer. 

Through  appropriate  programming  of  the  PET,  the  pressure  scanner  is 
stepped  to  the  next  pressure  line  to  repeat  the  measurement  process.  If 
probe  positioning  is  instead  required,  the  PET  activates  the  positioner 
by  closing  specific  channels  on  the  voltage  scanner.  The  collected  data 
are  processed  and  then  recorded  on  either  magnetic  cassette  tape  or  a 
disk  storage  system.  Data  may  also  be  viewed  on  the  CRT  of  the  PET  or 
printed  out  as  hard  copy. 

The  pressure  reference  system  (Ref.  2)  indicated  in  Fig.  5  provides 
high  precision  reference  pressures  for  on-line  calibration  of  the  trans¬ 
ducer  in  all  test  pressure  measurements.  Also,  through  the  fluid  switcli 
indicated,  different  pressures  may  be  supplied  at  various  steps  of  the 
pressure  scanner  as  reference  for  other  pressure  measurements.  Necessary 
temperature  measurements  are  also  made  using  thermocouples. 
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Figure  5.  Functional  diagram  of  data  acquisition  and 
experiment  control  system.  Task  I. 


To  illustrate  the  data  system,  a  program  designed  to  position  a 
five-hole  directional  pressure  probe  in  translation,  to  automatically 
null  the  probe  in  yaw,  and  to  record  probe  pressure  measurements  along 
with  associated  pressure  and  temperature  measurements  is  described  in 
the  flow  diagram  in  Fig.  6.  Raw  pressure  and  temperature  data  only  are 
collected  by  the  program;  further  reductions  of  data  to  obtain  total 
pressure,  Mach  number,  flow  pitch,  etc.,  based  on  the  five-hole  probe 
calibrations  (Ref.  3)  are  handled  in  other  programs  not  described  here. 
Pitot  tube  measurements  in  addition  to  the  five-hole  probe  measurements 
are  made  in  the  program  to  obtain  particular  reference  values  for  the 
flow  as  needed  in  data  reduction. 

The  program  is  made  up  of  five  major  parts  as  implied  in  the  flow 
diagram  by  the  five  branches  from  the  main  program.  These  parts  are: 

(1)  calculation  of  five-hole  probe  positions  for  translations,  (2)  pitot 
tube  measurements,  (3)  probe  translation,  (4)  automatic  yaw-nulling  of 
probe,  and  (5)  collection  of  probe  pressure  and  temperature  measurements. 
Operation  of  the  program  is  essentially  lef t-to-r ight  across  the  flow 
diagram,  with  positioning  of  the  probe  and  measurement  of  pressures 
(shown  toward  the  right-hand  portion  of  the  diagram)  carried  out  sequen¬ 
tially  over  all  the  previously  calculated  probe  positions.  On-line 
calibration  of  the  pressure  transducer  for  pressure  vs.  voltage  carried 
out  for  each  spin  of  the  pressure  scanner  provides  high  accuracy  in  all 
pressure  measurements. 

The  data  system  involving  the  PET  computer  as  just  outlined,  is  a 
self-supporting  system  entirely  controllable  through  appropriate  pro¬ 
gramming.  In  addition,  through  use  of  an  RS232  serial  adaptor  (see 
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again  Fig.  5)  the  PET  can  be  operated  as  a  remote  terminal  in  the  University 
Computational  Center  DATANET  system.  This  latter  system,  incorporating  a 
PDP  11/34  computer  and  high-capacity  disk  units,  is  a  university-wide 
high-speed  recirculating  digital  network  of  mini-  and  microcomputers  in 
experiment  control  and  data  acquisition.  A  file  management  system  called 
(UNIX)  (see  Ref.  4),  is  used  for  collection  of  experimental  data  and 
program  files,  and  assembled  object  code  can  be  directly  down-loaded  into 
the  remote  terminals  for  subsequent  execution.  Furthermore,  a  data  link 
exists  between  DATANET  and  the  University  mainframe  computer  which  pro¬ 
vides  an  even  further  expanded  data  reduction  and  computing  capability. 

1.4.  Presentation  and  Discussion  of  Experimental  Results 

Experimental  results  obtained  with  the  test  passage  are  presented 
in  Figs.  7  through  10.  These  plots,  based  on  preliminary  results  of  testing 
carried  out  to  date,  are  for  selected  traverse  planes  through  the 
passage.  Contour  plots  of  total  pressure  coefficient,  static  pressure 
coefficient,  and  Mach  number  are  shown  along  with  vector  plots  of  in¬ 
plane  (traverse  plane)  velocities.  In  traversing  the  five-hole  pressure 
probe,  21  points  were  measured  from  endwall  to  endwall  to  record  yaw  and 
pitch  angles  of  the  flow  and  impact  pressure.  This  procedure  was  carried 
out  for  8  or  9  slide  bar  positions  in  each  plane,  with  these  positions 
concentrated  mainly  toward  the  suction  sidewall.  Static  (or  total) 
pressures  obtained  with  the  sidewall  and  endwall  taps  were  combined  with 
the  traverse  data  to  complete  the  experimental  results  per  traverse  plane. 

A  summary  of  the  data  collection,  reduction  formulas,  coefficient  defi¬ 
nitions  and  coordinate  system  used  is  contained  in  Appendix  A.  In 
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generating  the  contour  plots  shown,  enrichment  of  the  measured  data  was 
employed  by  interpolating  at  two  intermediate  points  in  both  the  x  and  z 
coordinate  directions,  and  also  by  carrying  out  extrapolation  into  the 
corners  of  the  traverse  planes. 

Since  the  velocity  profile  generated  from  endwall  to  endwall  enter¬ 
ing  the  test  passage  was  symmetric  about  the  center  plane  of  the  passage, 
results  are  shown  only  for  the  lower  half  of  each  traverse  plane.  The 
entering  sidewall  boundary  layers  were  thin,  while  the  entering  endwall 
boundary  layers  were  relatively  thick,  producing  strong  vortical  motions 
in  the  main  stream  direction  throughout  the  passage.  The  movement  of 
low  energy  fluid  as  the  flow  proceeds  through  the  bend  and  discharge 
sections  can  be  traced  readily  in  terms  of  successive  plots  such  shown 
in  Figs.  7  through  10.  Also  flow  visualizations  were  made  by  producing 
ink  traces  on  the  lower  endwall  and  both  sidewalls  to  help  in  analyzing  the 
flow  and  to  guide  in  interpreting  measured  probe  data.  Such  flow 
visualization  produces  the  patterns  of  limiting  streamlines  on  the  passage 
walls  and  in  many  cases  may  or  not  indicate  flow  directions  interior  to 
the  passage  close  to  the  walls  because  of  possible  skewing  in  the 
boundary  layers  in  regions  of  high  shear.  These  patterns  showed,  as 
expected,  very  strong  cross-flow  on  the  endwalls  from  pressure  side  to 
suction  side  with  consequent  interaction  with  the  flow  on  the  suction 
sidewall  along  a  separation  line.  However,  no  evidence  was  seen  of 
three-dimensional  separation  of  the  entering  boundary  layer  itself  from 
the  endwall  which  has  been  found  to  predominate  the  flow  through  cascades 
(See  Refs.  5,  6). 
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The  plots  shown  in  Figs.  7  through  10  are  of  the  flow  in  selected 
traverse  planes  viewed  from  downstream,  with  the  suction  sidewall  on 
the  left  and  the  pressure  sidewall  on  the  right.  The  traverse  planes 
represented  are  the  inlet,  15-deg,  30-deg  and  90-deg  planes  through  the 
bend.  In  Fig.  7  the  development  of  total  pressure  loss  coefficient,  C 
can  be  seen.  The  higher  values  of  total  pressure  loss  corresponds  to 
higher  values  of  C  or  contour  label.  In  the  inlet  plane  (plane  1)  the 
Bernoulli  surfaces  in  the  entering  collateral  endwall  boundary  layer 
are  evident,  with  some  effect  of  the  bend  already  being  felt  in  growing 
of  the  pressure  side  boundary  layer.  Nothing  is  concluded  about  the  flow 
on  the  suction  sidewall  at  this  point,  since  the  limiting  position  of 
slide  bar  1  prevented  traversing  the  five-hole  probe  close  enough  to  the 
suction  side.  The  contours  of  C  in  traverse  planes  2  and  3  are  verv 

pt 

similar  to  those  in  plane  1,  showing  little  development  of  losses  through 
the  first  30  degrees  of  the  bend.  By  plane  7,  at  90  degrees  in  bend,  we 
can  see  an  increase  in  loss  level  along  with  thinning  of  the  endwall  layer 
in  combination  with  the  strong  interaction  of  the  endwall  and  sidewall 
flows  high  on  the  suction  side.  Also  in  plane  7,  the  pressure  sidewall 
layer  has  by  this  time  apparently  again  decreased  in  thickness. 

Figure  8  shows  contour  plots  of  static  pressure  coefficient,  in  which 
local  static  pressure,  pg,  is  referenced  to  the  pitot  static  and  dynamic 
pressures  pgQ  and  qQ,  respectively,  for  the  entering  velocity  profile. 
These  contours,  running  vertically  from  endwall  to  endwall,  correspond 
to  those  for  potential  flow  over  the  great  majority  of  the  passage  for 
planes  1  through  3,  but  with  increasing  pressure  gradient  across  the 
passage  from  suction  side  to  pressure  side.  The  center  of  low  pressure. 
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in  plane  7,  lias  moved  off  the  section  sidewall  out  into  Lhe  flow  along 
with  formation  ot  tin'  passage  vortex  near  tlie  suction  side.  The  overall 
gradient  of  pressure  from  sidewall  to  sidewall  has  now  decreased  in 
response  to  the  action  of  the  losses  in  the  flow.  Also,  an  extraneous 
vortex  filament  is  evident  from  the  local  low  pressure  center  near  tin- 
pressure  sidewall.  Further  testing  of  the  flow  is  needed  to  explain  this. 

In  Fig.  9  we  see  the  in-plane  velocity  component  in  the  four  traverse 
planes.  For  the  inlet  plane,  although  these  components  are  small,  on  the 
order  of  0  to  5  or  6  ft/s,  there  is  indication  of  downward  flow  on  the 
pressure  side  and  across  the  endwall  toward  the  suction  side,  with  possible 
starting  of  the  passage  vortex  in  the  suction  side  corner  at  the  endwall . 

In  planes  2  and  3  the  probe  traverse  closest  to  the  suction  sidewall  is 
still  outside  the  boundary  there  as  indicated  by  the  in-plane  velocity 
vectors  pointing  essentially  outward  from  the  wall.  Also  in  plane  3 
there  is  a  slight  indication  of  flow  across  the  centerline  from  the  upper 
half  of  the  passage.  For  planes  2,  3  and  7  the  strong  secondary  flow 
formation  in  development  of  the  passage  vortex  is  seen. 

In  Fig.  10  of  Mach  number  contours,  planes  1,  2,  and  3  show  limited 
secondary  flow  action,  except  for  the  spanwise  velocity  gradients  caused 
by  the  passage  curvature.  In  plane  7,  however,  the  low  energy  region  in 
the  suction-side  corner  is  being  wrapped  up  into  the  passage  vortex  that 
forces  itself  down  into  the  corner.  Again  we  see  evidence  of  the  extraneous 
vortex  near  the  pressure-side  corner.  Plane  7  also  shows  an  overall 
increase  in  average  velocity  which  is  caused  by  the  convergence  of  the 
passage  cross  section. 
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TASK  11:  KVALIA1  ION  AND  MANACKMKNT  OK  UNSTKADY  FLOW  F.FFL(TS  IN 
AXIAL- FLOW  COMPKKSSOkS 

In  order  to  further  our  knowledge  and  experience  involving  unsteady 
flow  effects  in  axial-flow  compressors  acquired  under  AI'OSR  grant  76-2916  a 
detailed  review  of  related  literature  was  started  in  the  Kail  of  1978  and 
continues.  Also,  more  experiments  were  completed  on  the  three-stage  research 
compressor  used  previously  (see  Figure  11).  The  conclusions  reached  by  the 
beginning  of  1979  are  documented  in  Reference  7.  Subsequent  progress  in 
Task  IT  achieved  under  AKOSR  Contract  F49620- 79-C-0002  is  summarized  in  Lhe 
following  paragraphs. 

A  note  about  the  data  that  will  be  described  below  may  be  helpful. 
Time-average  results  were  acquired  with  slow  responding  instruments  and 
represent  the  average  value,  over  time,  of  a  continuous] v  measured  quantity. 
Periodic-average  results  involve  Lhe  average  value  of  a  periodical ly  (e.g. 
once  per  rev)  sampled  quantity.  Blade-to-bl ade-average  results  are  the 
spatially  integrated  values  of  blade-to-blade  profiles  of  either  time-average 
o£  periodic-average  quantities.  Time-average  results  are  independent  of 
rotor  sampling  position  while  periodic-average  results  are  generally  not. 

To  ascertain  whether  or  not  varying  the  relative  circumferential 
positions  of  the  blades  (same  number  in  each  row)  of  the  inlet  guide  vane 
(1GV)  and  stator  rows  of  the  three-stage  research  compressor  would  result  in 
measurable  time-average  stator  performance  (deviation  angle  and  total-head 
loss)  changes,  the  compressor  was  operated  at  2400  rpm  and  local  time-average 
total-head  measurements  were  acquired  for  two  different  circumferential 
position  settings  of  the  ICV  and  first  stator  rows.  At  a  lower  rotor  speed 
(1400  rpm),  similar  tests  indicated  an  appreciable  difference  in  inleL  noise 
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Figure  11.  Research  compressor  with  probe  measurement  stations 
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(about  6dB  in  the  octavo  bands  nearest  Lin*  blade  passing  frequency)  for 
different  relative  circumferential  positions  of  the  iGV  and  stator  rows. 


However,  time-average  total -head  loss  and  deviation  angle  data  were 
inconclusive  about  stator  performance  differences.  At  the  higher  rotor 
speed,  minimum  and  maximum  noise  circumferential  positions  of  the  IGV  and 
first  stator  rows  involved  a  difference  in  noise-sound  pressure  level  of 
about  10  dB  in  the  octave  band  including  the  blade  passing  frequency. 

Contour  graphs  of  the  time-average  total-head  measurements  behind  the  first 
stator  are  displayed  in  Figure  12.  Blade-to-blade-average  values  of  the 
information  in  Figure  12  at  various  radii  are  compared  in  Figure  13. 
Representative  time-average  yaw  angle  data  are  shown  in  Figure  14.  The  blade 
and  measurement  position  parameters  are  explained  in  Figure  15.  Although  the 
annulus  cross-section  distribution  of  stator  exit  local  time-average  total 
head  varies  appreciably  with  change  of  relative  circumferential  positioning 
of  IGV  and  stator  blades  the  blade-to-blade-average  total-head  values  do  not. 
Also,  the  stator  exit  free-stream  angles  do  not  seem  to  depend  significantly 
on  circumferential  positioning  of  the  IGV  and  stator  blades.  Thus  it 
appears  as  if  the  circumferential  placement  of  IGV  and  first  stator  blade 
rows  relative  to  each  other  did  not  result  in  measurable  stator  performance 
differences  under  the  present  test  conditions. 

A  fast-response  total-pressure  probe  was  used  with  a  periodically 
sampling  and  averaging  measurement  system  to  observe  the  details  of  unsteady 
total-pressure  variations  in  the  first  stage  discussed  above.  The  compressrr 
was  operated  at  the  minimum-noise,  1400  rpm  condition  used  in  earlier 
studies.  Details  of  the  research  are  documented  in  Reference  8.  In  general, 
periodic-average  total -pressures  varied  appreciably  with  rotor  sampling 
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Figure  13.  Blade-to-blade -average  first  stator  exit  flow  total-head 

variation  with  stator/lGV  circumf erential  positioning  (2400 


-AVERAGE  FLOW  ANGLE 


CIRCUMFERENTIAL  POSITION  OF  MEASUREMENT  Y/Ss 
(a)  YOiqv/Ss  =  0;  Y01$/Ss  =  0.159 


t/D 

UJ 


-J  sc 
U-  o 


CD  o 


UJ 

5*  O 
»  — J 


CIRCUMFERENTIAL  POSITION  OF  MEASUREMENT  Y/S$ 


(b)  Y0IGV/S$  =  0;  Y0ls/S$  =  0.719 


Figure  14.  Tirae-average  first  stator  exit  flow  angle  variation  with 

stator/IGV  circumferential  positioning  at  50  percent  passage 
height  from  the  hub  (2400  rpm) . 
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Figure  15.  Blade  cascade  showing  relative  positions  of  blades  for  several  rotor  sampling 
positions.  (SR  is  rotor  blade  spacing;  Ss  is  stator  blade  spacing;  Y  is 
measurement  circumferential  location;  Y0R  is  reference  rotor  blade  circum¬ 
ferential  location;  YO-^g  is  reference  first  stator  blade  circumferential 
location. ) 


position  downstream  of  the  first  and  second  rotors  buL  very  little  downstream 
of  the  first  stator.  Some  of  the  daLa  acquired  at  mid-span  that  demonstrate 
this  observation  are  shown  in  figures  16  and  17.  Figures  16(a),  (b)  and  (c) 
involve  comparisons  of  local  periodic-average  and  time-average  total  pressures. 
Figure  17  involves  a  comparison  of  blade-to-blade  averages  of  the  periodic- 
average  and  time-average  total-head  profiles  shown  in  Figure  16. 

Similar  behavior  was  noted  at  other  span  locations.  The  Lrends  in  total  - 
pressure  unsteadiness  could  be  explained  in  terms  of  the  blade-wake  chopping, 
transport  and  interaction  kinematics  model  proposed  in  Reference  9.  Since 
total-pressure  values  are  indicative  of  fluid  particle  energy  addition 
through  work  and  energy  loss  through  friction.  Table  1  could  be  organized 
as  shown  to  help  explain  what  was  observed  downstream  of  rotors  1  and  2.  The 
periodic-average  total-head  profile  shape  and  magnitude  downstream  of  the 
rotors  seemed  to  reflect  the  different  kinds  of  particles  of  fLuid  identified 
in  Table  A  chat  were  in  the  measurement  window  (between  the  +'s  in  Figure  15) 
for  any  rotor  sampling  position.  Periodically  varying  spatial  distributions 
and  proportions  of  the  different  kinds  of  fluid  particles  at  the  measurement 
window  resulted  in  differing  total-head  profiles.  This  kind  of  periodic 
fluctuation  was  not  as  prominent  behind  the  stator,  largely  because  all 
fluid  particles  lost  energy  only  in  the  stator  row. 
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(a)  first  rotor  exit,  midspan 


Table  1.  Behavior  of  fluid  particles  moving  from  the  rotor  inlet  to 
the  measurement  station  downstream  of  the  rotor 


Type  of 
fluid 
particle 

Total  pressure 
at  the 
rotor  inlet 

Energy 

added 

Total  pressure 
at  the  rotor  exit 
Losses  measurement  station 

Freestream 

A 

B 

C 

D 

Chopped 

IGV  wake 

<A 

>B 

>C 

>D 

Noninteracted 
rotor  wake 

A 

»B 

>C(50%-90%  PHH) 

>D 

»C(30%) 

D 

»>C(10%) 

<D 

Interacted 
rotor  wake 

<A 

»B 

»>c 

«D 

The  freestream  values  are  set 
values  are  estimated  in  terms 
than  these  reference  values. 

arbitrarily  at 
of  being  less 

levels  A,  B,  C  and  D.  Other 
than,  equal  to,  or  greater 

3.  TASK  III:  DEVELOPMENT  OF  A  UKSICN- SYSTEM  ORIENTED  DEVIATION  ANTED 
PREDICTION  EQUATION  FOR  ADVANCED  COMPRESSOR  AND  FAN  CON FI  CURAT  I ONS 

Earlier  analytical  (USAF/AFAPE  Contract  F33615-76-C-2090)  [Ref.  10 | 
and  computational  (NASA  (.'rant  NsG-3033)  [Ref.  11]  work  suggested  a  correla¬ 
tion  of  compressor  and  fan  cascade-flow  turning  angles  with  blade-profile 
boundary  layer  development  characteristics.  The  results  of  these  investi¬ 
gations  were  integrated  and  extended  (USAF/AFOSR  Grant  78-3609)  to  cover 
more  challenging  test  cases  including  more  of  the  variables  existing  in 
real  compressor  flows.  Task  III  of  the  current  contract  was  directed 
toward  the  development  of  a  design-system  quality  deviation  angle  prediction 
equation  for  advanced  compressor  and  fan  cascade  configurations. 

It  was  found  in  Reference  10  that  there  were  few  cascade  data  sets 
available  with  adequate  information  to  give  insight  into  Lhe  nature  of  the 
interaction  of  the  boundary  layer  characteristics  and  the  blade  cascade 
deviation  angle.  Therefore,  the  computation  method  of  References  11  and  12 
was  used  in  the  current  investigation  to  predict  the  blade-to-blade  passage 
boundary  layer  development  characteristics  and  the  outlet  fluid  flow  angles 
for  typical  cascade  geometries.  Test  configurations  which  involved  many  of 
the  variables  existing  in  real  compressor  flows  were  chosen.  Each  test 
case  included  some  experimental  data  which  was  used  to  verify  Lhe  validity 
of  the  computational  model. 

Several  changes  were  made  in  the  original  computer  program  to  better 
model  the  additional  test  configurations.  A  new  differentiation  scheme  was 
used  to  determine  the  surface  pressures  near  the  trailing  edge  of  the  blade. 
This  differentiation  method  provided  smoother  and  more  consistent  surface 
pressure  distributions  than  did  the  previous  spline  f it t ing *method .  The 
better  estimation  of  surface  pressures  also  improved  calculation  convergence 
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Tin*  program  was  also  modified  to  modol  small  regions  of  supersonic 
flow.  The  modification,  suggested  and  outlined  by  Jerry  Wood  of  the 
NASA-l.ewis  Research  Center,  has  proven  to  he  useful,  though  when  compared 
in  experimental  data  iL  underestimates  the  peak  Mach  number.  The  total - 
pressure  loss  model  of  Reference  I!  and  12  was  also  modified  to  account 
for  changes  in  the  Lhickness  and  the  radial  location  of  the  streamsheet 
downstream  of  the  trailing  edge. 

A  total  of  eight  Lest  cascade  configurations  involving  a  wide  range 
of  variables  were  selected.  The  variables  included  incidence  angle,  inlet 
Mach  number,  inlet  Reynolds  number,  and  axial-veloc i l v  density  ratio  (AVDR) 
Table  2  shows  the  test  configurations  and  the  significant  variables  for 
which  computations  were  made. 

The  double  circular  arc  cascade  [Ref.  13]  provided  insight  into  the 
effect  of  the  axial  velocity-density  product  ratio  (AVDR)  on  the  flow 
field.  Figures  18  and  19  show  the  large  discrepancy  which  can  exist  between 
the  surface  pressures  if  the  distribution  of  AVDR  is  not  correctly  modeled. 
Additional  insight  into  the  effect  and  distribution  of  AVDR  was  gained 
through  discussions  with  H.  Starken  of  the  DFVLR.  The  importance  of 
determining  whether  the  AVDR  reported  in  experiments  corresponds  to  a 
ratio  taken  between  exit  and  inlet  measuring  stations  or  between  trailing 
and  leading  edge  averages  was  verified  by  computation.  This  ratio  has  a 
particularly  significant  effect  on  pressure  distribution,  turning  and  loss 
for  high-diffusion-ratio  cascades. 

Variation  of  incidence  for  the  10A30/27 . 6tt4  5  cascade  [Ref.  14]  pro¬ 
duced  computed  fluid  turning  which  was  similar  to  that  experimentally 


measured.  The  experimental  measurements  also  included  the  distribution  of 
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Figure  19.  Calculated  Surface  Pressures  Compared  with  Experimental 
Data  for  the  DCA  <)>  =  48°  Blade  Cascade. 

Mi  =  0.64  Rec  =  1.26  x  10&  ^  =40°  i  =  1° 

D  =  0.244  y  =  15°  0  =  1.629  c  =  1.24 


47 


velocity  throughout  ihe  flow  field.  Very  s.-u  isfanorv  flow  fi  Id  eompar  i  - 
sons  between  experiment  and  computation  were  reported  at  an  inlet  Mach 
number  of  0.41  in  Reference  12.  At  an  inlet  Mach  number  of  0.79,  computed 
maximum  flow  field  velocities  were  lower,  but  si  i  I  1  favorable,  vlien 
compared  wiLli  the  experimental  measurements. 

The  computed  flow  through  the  65-(12)10  cascade  [Rei.  I’]  at  low 
Reynolds  number  displays  the  effect  of  the  boundary  layer  on  direction  of 
fluid  flow  (Figure  20).  Computed  turning  at  high  Mach  numbers  did  not:  compare 
favorably  with  Lite  measurements,  probablv  because  axial -velocil v  density 
ratios  other  than  1.0  occurred  in  Lite  experimental  data. 

Tiie  computed  high-speed  flow  through  Lite  65-  (1 2A.,  1^,  )  10  cascade 
[Ref.  15]  compared  well  with  experimental  measurements.  At  higher  Mach 
numbers,  there  was  again  an  underestimation  of  the  peak  Mach  number. 

The  0NERA  115  Stator  cascade  [Ref.  16)  was  the  only  configuration 
for  which  blade-surface  boundary  layer  measurements  were  available.  The 
comparison  of  the  calculated  boundary  displacement  and  momentum  thickness 
distribution  with  the  experimental  measurements  was  very  good  (Figure  21). 

The  lack  of  cascade  configuration  boundary  layer  development  data  is  an 
important  weakness  in  support  for  computation  method  development  and  for 
turning  angle  prediction  equation  improvement. 

The  Prescribed  Velocity  Distribution  (PVD)  cascade  [Ref.  17]  provided 
a  series  of  test  cases  at  low  Reynolds  number  for  which  ihe  AVDK  was  varied 
at  several  incidence  angles.  These  data  and  Lite  corresponding  computation 
points  showed  the  change  in  the  boundary  laver  characteristics  with  the 
variation  of  incidence  and  axi al -velocity  density  ratio. 

The  1,0504  test  cascade  has  considerable  potential  for  useful  results. 

The  cascade  geometry  represents  a  cascade  plane  projection  of  a  section 
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Figure  20.  Calculated  Flow  Angle  at  the  Edge  of  the  Suction  and  Pressure 
Surface  Boundary  Layers  for  the  NACA  65-(12)10  Blade  Cascade. 

*1  -  0.06,  Rec  -  1.47  *  10^,  ^  -  57.7°,  a  =  12°, 

D  «  0.42,  y  -  45.7°,  a  ■=  1.0,  c  ■=  124  mm,  AVDR  =  1.0 


from  the  rotor  of  Reference  18.  As  a  part  of  the  test  case  program  of 
AGARi)  Working  Group  12  on  Through-Flow  Calculation  in  Turbomachines,  rotor 
laser  two-focus  measurements  were  released  along  with  linear  cascade 
performance  measurements.  The  cascade  information  included  high  Mach  number 
turning  and  loss  for  various  inlet  flow  angle,  with  one  data  point  at  the 
choked- flow  condition.  Flow  field,  turning  and  loss  computation  runs  were 
completed  at  each  of  the  three  data  points.  The  results  have  not  been 
evaluated  at  this  time,  but  will  be  studied  in  relation  to  other  AGARD 
WO.12  calculations  made  in  France  and  Belgium. 

A  NACA  65-(18)15  cascade  [Ref.  19]  was  experimentally  tested  both  as 
a  radial  and  an  axial  cascade.  Computations  were  made  for  both  radial  and 
axial  cascade  arrangements  to  show  the  effect  of  the  radial-  and  axial 
velocity  density  ratio  on  the  cascade.  The  results  were  reported  in 
Reference  19. 
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4.  TASK  TV:  EXPERIMENTAL  STUDY  OF  MERITS  OF  SPECIFIC  FLOW  PATH  GEOMETRY 
CHANCES  IN  CONTROLLING  SECONDARY  FLOWS  IN  AXIAL-FLOW  COMPRESSORS 

The  suitability  of  using  the  existing  Iowa  State  research  compressor 
blading  (used,  for  example,  in  Task  11)  for  a  "baseline"  configuration  for 
performance  comparison  with  a  "modified"  configuration  was  examined  at  the 
beginning  of  Lhe  contract  period.  It  was  decided  that  a  new  set  of  baseline 
blades  should  be  designed  and  built.  Higher  blade-c.hord  Reynolds  number, 
more  suitable  blade  material,  conventional  blade  section  shapes,  a  favorable 
ratio  of  number  of  rotor  blades  to  number  of  stator  blades,  and  elimination 
of  inlet  guide  vanes  were  some  of  the  improvements  expected. 

The  new  baseline  blades  were  designed  with  the  aid  of  computer  codes 
described  in  References  20  and  21.  The  preliminary  design  code  [Ref.  20] 
is  based  on  simple  radial  equilibrium  of  the  duct  flow  between  blade  rows. 

A  constant  spanwise  distribution  of  blade-element  efficiency  was  involved 
for  each  blade  row.  For  the  preliminary  design  phase,  a  free-vortex  swirl 
velocity  distribution  behind  the  rotor  was  used.  Values  of  input  parameters, 
including  rotor  tip  and  stator  hub  diffusion  factors,  rotor  and  stator  tip 
axia 1 -veloc iLy  ratios,  blockage  factors  and  rotor  and  stage  polytropic 
efficiencies,  were  varied  until  reasonable  input  and  output  were  achieved. 

The  final  design  code  [Ref.  21]  is  based  on  a  streamline  curvature  type 
solution  of  Lhe  two-dimensional  axisvmmetric  approximation  of  actual 
axial-flow  compressor  flow.  It  involves  duct-flow  stations  outside  blade 
rows.  VeLocity  diagrams,  however,  are  determined  at  blade  leading  and 
trailing  edges.  The  code  produces  detailed  aerodynamic,  as  well  as  blade 
geometry  data. 

The  new  baseline  configuration  consists  of  two  identical  stages  as 
illustrated  in  the  meridional  plane  sketch  of  Figure  22.  The  flow  path 
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Figure  22.  Meridional  plane  view  of  new  blading. 


dimensions  involved,  as  well  as  the  rotor  drive  system  and  flow  rate  control, 
are  the  same  as  before  as  described,  for  example,  in  Ref.  9.  A  summary  of 
pertinent  design  data  is  given  in  Table  3.  Representative  blade  section 
profiles  are  shown  in  Figures  23  and  24.  Blade-element  losses  were  approxi¬ 
mated  with  values  from  Figure  203  of  Reference  22.  Blade-element  incidence 
and  deviation  angles  were  estimated  using  the  method  proposed  in  Chapter  VI 
of  Reference  22.  Approximate  blockage  factors  used  during  the  preliminary 
design  phase  were  specified.  Constant  energy  addition  over  the  rotor  span 
was  selected. 


Table  3.  Summary  of  Baseline  Compressor  Design  Data 


Rotor  speed 

2400  rpm 

Flow  rate 

5.25  lb  /s  (2.38  kg/s) 

Pressure  ratio 

1.019 

Number  of  blades 

Rotor 

21 

Stator 

30 

Blade  material 

fiberglass  (steel  trunnion) 

Blade  chord 

2.39  in.  (6.07  cm) 

Blade  profile 

double  circular  arc 

Flow  path 

Hub  radius 

5.60  in.  (14.22  cm) 

Tip  radius 

8.00  in.  (20.32  cm) 

A  modified  stator  blade  with  forward  symmetrical  sweep  of  the  blade 
leading  edge  in  the  endwall  regions  was  designed.  Air  Force  Aero  Propulsion 
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Representative  baseline  configuration  stator  blade  sections. 
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Laboratory  experience  suggests  tliat  such  a  modification  may  lead  to  improved 
stator  performance.  The  previously  mentioned  computer  code  [Ref.  2 1 ]  was 
used  to  aerodynamicallv  design  this  blade  for  use  with  the  baseline  rotor. 
Except  for  nonconstant  chord  distribution,  the  baseline  design  options  were 
used.  Some  representative  blade  section  profiles  are  shown  in  Figure  >. 
Because  of  the  unusual  nature  of  this  blade,  an  independent  analysis  was 
made  using  the  AFAPL  UD0300  axial-flow  compressor  design  code  [Rei>.  2  3  and 
24}.  Like  the  program  in  Reference  21,  UD0300  employs  the  streamline 
curvature  method  to  determine  the  axisymmetric  approximation  of  the  actual 
compressor  flow.  Unlike  the  program  in  Reference  21,  UD0300  allows  place¬ 
ment  of  calculation  stations  within  blade  rows.  A  comparison  of  incidence 
angles  predicted  by  both  codes  for  the  modified  blade  conf igurat ion  designed 
with  the  code  of  Ref.  21  is  shown  in  Figure  26.  Since  the  incidence  angles 
were  similar,  it  was  decided  that  the  modified  blade  configuration  designed 
with  the  NASA  code  was  suitable  for  use. 

Performance  testing  of  the  baseline,  compressor  is  scheduled  to  begin  in 
late  1980  and  will  consist  mainly  of  time-average  and  periodic-average 
total-pressure  and  flow  velocity  measurements.  Design-point  operation  flow 
field  properties  will  be  ascertained  in  detail.  Some  off-design  data,  such 
as  the  rotating  stall  limit,  will  be  acquired.  The  baseline  configuration 
blades  have  been  fabricated  and  have  been  mounted  on  the  rotor  hub  and  stator 
rings.  Blade  tip  clearances  are  being  machined  and  the  baseline  compressor 
testing  should  commence  in  December  1980.  Eyelash  gage  checks  of  the  blade 
profiles  indicate  very  good  correspondence  with  the  design  coordinates. 

Testing  of  the  modified  compressor  is  anticipated  to  begin  in  1981. 
Fabrication  of  the  modified  stator  blades  will  begin  soon;  the  contract 
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Figure  26.  Comparison  of  modified  stator  leading  edge  incidence  angle  predictions  calcu 
using  different  compressor  design  computer  codes. 


involved  has  been  awarded.  Time-average  and  periodic-average  total-pressur 
and  t low  velocity  measurements  for  design-point  and  off-design  operation  wi 
be  acquired  and  compared  with  baseline  results. 
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5.  TASK  V:  DEFINITION  OF  EXPERIMENTAL  PROGRAMS  AND  FACILITIES  APPROPRIATE 
FOR  UNIVERSITY  TURBOMACHINERY  RESEARCH  PROGRAMS 

A  very  wide  range  of  research  and  development  activity  is  currently 
carried  on  in  the  field  of  turbomachine  aerodynamics.  This  research  and 
development  takes  place  in  all  highly-industrialized  countries  of  the  world, 
and  in  a  few  countries  that  are  less  developed.  Some  characteristics  of 
turbomachinery  research,  as  distinguished  from  development,  should  be 
recognized : 

1.  In  the  United  States,  for  many  years,  the  position  was  taken  in 
important  elements  of  the  aircraft  propulsion  and  industrial 
turbomachinery  industries  that,  because  of  the  "proprietary" 
nature  of  the  information  generated  in  turbomachinery-relaied 
experiments,  turbomachine  aerodynamic  research  was  industrial 
business  and  not  appropriate  for  university  programs.  To  a  lesser 
degree,  and  for  slightly  different  reasons,  some  government 
research  funding  agencies  took  the  same  point  of  view.  To  some 
extent,  this  attitude  exists  today. 

2.  In  the  remainder  of  the  world,  a  different  view  exists.  Every 
country  which  is  technically  competitive  with  the  United  States 
has  from  two  to  several  active  university  research  programs  in 
turbomachinery  fluid  dynamics,  usually  in  a  formally  organized  and 
supported  institute  or  laboratory  with  a  Professor  as  director, 
and  a  technical  staff. 

3.  The  facilities  available  to  university  turbomachinery  research 
groups  in  the  United  States  are  modest,  marginal  or  inadequate. 

A  number  of  "apparent"  research  programs  are  largely  based  on  work 
generated  through  industrial  consulting  or  on  data  supplied  from 
government  laboratory  in-house  experiments. 
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Turbonuuh  i  nor  v  exper  i men La  1  facilities  in  European  and  Japanese 


universities  are  custom-designed  and  highly  appropriate  for 
continuing  research.  Most  are  comparable  in  quality  to  industrial 
and  government  research  laboratory  facilities  in  the  United  States 
Some  turbomach inary  aerodynamic  research  should  not  be  carried  on  in 
universities.  Prohibitive  facility  operating  cost  and  safety  consideration; 
might  be  typical  reasons  for  disqualification  of  some  classes  of  research. 
However,  there  are  several  categories  of  research  projects  which  can  be 
considered  as  well-matched  to  university  characteristics,  including  faculty 
qualifications,  student  capability  and  available  supporL  services.  ihree 
are  outlined  in  Sections  1.1  to  5.3. 


5.1.  Experiments  in  Low-Speed  Rotating  Test  Facilities 

Performance  significant  aerodynamic  phenomena  in  fan,  compressor  and 
turbine  flow  path  components  can  be  effectively  studied  in  low-rot  at ional - 
speed  test  units.  At  moderate  to  large  scale,  useful  detailed  measurements 
can  be  made  at  moderate  to  low  cost. 

Secondary  flows  and  geometrical  means  for  modifying  or  controlling 
secondary  flows  remain  an  unresolved  problem  area  in  axial  and  cenLrifug.il 
compressors  and  turbines.  The  careful  definition  of  three-dimensional 
effects  by  measurements  is  necessary  to  permit  design  improvements  and  also 
to  produce  test  or  "benchmark"  data  for  development  of  computational  mo i hod 
Present  availability  of  such  data  is  totally  inadequate  and  new  ideas  an 
genuinely  needed. 

The  widespread  significance  of  the  essential  unsteadiness  in  turbo- 
machinery  flows  has  become  Increasingly  evident.  Although  useful  work  has 
been  done  in  university  Lurbomachinerv  research  facilities  in  exposing 
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important  feaLures  of  unsteady  flows  in  compressors  and  turbines,  a  great 
deal  of  quantitative  evaluation  remains  to  be  done.  This  type  of  investiga¬ 
tion  is  being  carried  on  in  aircraft  engine  industry  laboratories,  but 
additional  large-scale,  low-speed  facilities  could  be  utilized  without 
overcrowding  the  field.  The  patient  and  original  work  of  advanced  students 
would  be  a  valuable  asset  in  both  experimentation  and  data  analysis. 

Nearly  all  currcnL  low-speed  compressor  and  turbine  test  facilities 
are  single-stage  configurations.  Both  the  unsteady  and  secondary  flow 
effects  should  also  be  studied  in  multistage  environments.  A  multistage 
test  unit  could  be  efficiently  and  economically  operated  for  research  purposes 
in  a  university. 

5.2.  Experiments  in  Linear  Cascade  Facilities 

The  computational  study  carried  out  under  Task  III  of  Lhe  present 
program  demonstrated  clearly  that  there  is  insufficient  information  available 
on  the  flow  fields  found  in  linear  airfoil  cascades.  To  permit  development 
of  more  satisfactory  turning  and  loss  prediction  systems  for  compressor  and 
turbine  cascade  flows,  a  carefully  laid  out  sequence  of  test  case  data 
points  should  be  made  available.  There  are  only  a  few  linear  cascade 
facilities  now  operational  in  the  United  States.  Only  one  approaching 
test  case  quality  exists  in  a  university.  This  contrasts  greatly  with  the 
situation  in  Europe  where  graduate  student  research  in  cascade  facilities 
has  generated  a  substantial  fraction  of  the  cascade  flow  field  information 
available  in  the  world  today. 

5.  1.  Experiments  Evaluating  and  Introducing  New  or  Improved  Measuring  Techniques 
A  Lliird  class  of  research  suited  to  university  programs  is  the  develop¬ 
ment  of  non-L rad i t ional  instrumentation  and  measurements  systems  for 

6  1 


turbomachinery  flows.  Detailed  flow  field  mapping,  and  measurements  in 
unsteady  turbomachine  flows  are  manpower-intensive  research  activities. 
Universities  now  engaged  in  turbomachinery  experimentation  have  been 
especially  successful  in  the  introduction  of  new  measurement  techniques. 

The  scope  of  typical  projects  in  measurement  methods  is  such  that  individual 
problems  can  be  isolated  and  attacked  readily  by  facul ty-studenL  groups. 
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SECTION  TV 
PROGRAM  PERSONNEL 

Three  principal  investigators  share  responsibility  for  the  current 
program. 

-  George  K.  Serovy,  Anson  Marston  Distinguished  Professor  in 

Engineering,  Tasks  III,  IV,  V 

-  Patrick  Kavanagh,  Professor  of  Mechanical  Engineering,  Task  I 

-  Theodore  H.  Okiishi,  Professor  of  Mechanical  Engineering,  Tasks  II, 

IV 

Four  graduate-level  engineers  were  also  associated  with  the  work. 

-  John  A.  McAndrew,  Graduate  Research  Assistant,  Task  I 

-  Elmer  C.  Hansen,  Post-Doctoral  Research  Fellow  (to  1  October  1979), 

Tasks  III,  IV 

-  William  C.  Zierke,  Graduate  Research  Assistant,  Task  II 

-  Michael  D.  Hathaway,  Graduate  Research  Assistant,  Task  IV 
During  the  research,  several  undergraduate  students  in  Mechanical 

Engineering  made  useful  contributions.  This  was  an  effective  mechanism  for 
introducing  undergraduates  to  the  research  activity  in  turbomachinery. 

-  Mark  T.  Diefenthaler ,  Undergraduate  Assistant,  Task  I 

-  Allen  D.  Dvorak,  Undergraduate  Assistant  (to  1  June  1980),  Task  I 

-  Douglas  A.  Hottman,  Undergraduate  Assistant,  Task  I 

-  J.  Scott  Meline,  Undergraduate  Assistant  (to  1  June  1979),  Task  I 

-  Betsy  D.  Morgan,  Undergraduate  Assistant  (to  1  June  1979),  Task  1 

-  Brian  H.  Pigg,  Undergraduate  Assistant  (to  1  June  1979),  Task  I 

-  Paul  G.  Smith,  Undergraduate  Assistant,  Task  I 


68 


r 


SECTION  V 

INTERACTION  WITH  FEDERAL  AGENCIES  AND  INDUSTRY 
The  turbomachinery  research  program  at  Iowa  State  University  has 
focused  on  projects  which  make  a  contribution  to  the  development  of  design 
systems  for  advanced  compressors,  fans,  and  turbines  for  air-breathing  air¬ 
craft  propulsion  systems.  The  current  contract  has  not  changed  this  focus 
and  has  involved  numerous  direct  contacts  with  outside  organizations. 

Task  II  is  a  continuation  of  turbomachine  unsteady  flow  measurement 
work  initiated  under  AFOSR  Grant  76-2916.  Related  interactions  involving 
formal  written  discussions  and  informal  conversations  about  measurement 
techniques  and  data  arc  summarized  below. 

Organization  and  Nature  of  Contact 
von  Karman  Institute  for  Fluid  Dynamics, 

Rhode-St.-Genese,  Belgium;  formal  lectures, 
lab  visit,  other  technical  discussions 

Institut  fur  Strahlantriebe  und 
Turboarbeitsmaschinen,  Technische 
Hochschule,  Aachen,  West  Germany;  lab 
visit  and  technical  discussions,  review 
of  technical  paper 

Turbopropulsion  Laboratory,  Department  of 
Aeronautics,  Naval  Post-graduate  School, 

Monterey,  California;  formal  lecture,  lab 
visit,  and  technical  discussions 
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Individual  Contacts 
F.A.E.  Breugelmans 


H.  E.  Gallus 


R.  P.  Shreeve 


NASA-Lewis  Research  Center,  Cleveland 


C.  Ball 


Ohio;  lab  visits  and  technical 

M.  Hartmann 

discussions 

A.  .1.  Strazisar 

The  Trane  Company,  LaCrosse,  Wisconsin; 

telecom  technical  discussions 

D.  B.  Schmidt 

United  Technologies  Research  Center, 

East  Hartford,  Connecticut;  telecom 

technical  discussions 

J.  H.  Wagner 

Department  of  Civil  and  Mechanical 

Engineering,  University  o'  Tasmania, 

Hobart,  Tasmania;  technical  discussions 

G.  J.  Walker 

Department  of  Aerospace  Engineering, 

B.  Lakshminaravana 

Pennsylvania  State  University,  University 
Park,  Pennsylvania;  technical  exchange  of 
data,  formal  discussion  of  research 

University  of  Tennessee  Space  Institute,  M.  Kurosaka 

Tullahoma,  Tennessee;  formal  lecture  J.  E.  Caruthers 

Task  III  has  involved  use  of  experimental  test  cases  selected  to  demon¬ 
strate  capabilities  and  limitations  of  the  cascade  flow  field  computation 
system.  It  has  also  involved  continued  cooperation  with  the  initial  sponsors 
at  the  NASA-Lewis  Research  Center.  Mr.  M.  J.  Hartmann,  Chief,  Fluid  System 
Component  Division,  Lewis  Research  Center,  has  approved  use  of  the  NASA 
computation  facilities  for  test  case  examples. 
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Organization  and  Nature  of  Contact 


Individual  Contacts 


J.  R.  Switzer 

Dr.  G.  Winterfield 
Dr.  H.  Starken 
Dr.  J.  Renken 
Dr.  H.  B.  Weyer 

J.  Fabri 
G.  Meauze 


M.  J.  Hartmann 
P.  M.  Sockol 
D.  M.  Sandercock 
J.  F.  Schmidt 
J .  Wood 

Task  IV  depends  on  substantial  cooperation  with  USAF/AFAPL  and  industry. 

A  list  of  the  visits  and  telephone  contacts  for  Task  IV  follows. 

Organization  and  Nature  of  Con tac t  Individual  Contacts 

NASA- Lew is  Research  Center,  Cleveland,  D.  M.  Sandercock 

Ohio;  exchange  of  computer  programs,  J.  E.  Crouse 

review  of  designs 

I 


Air  Research,  Phoenix  Division,  Phoenix, 
Arizona;  test  cases 

Institut  fur  Antriebstechnik,  Deutsche 
Forschungs-  und  Versuchsanstal t  fur  Luft- 
und  Raumfahrt,  Koln,  West  Germany;  test 
cases,  review  of  results 

Direction  de  1 ' Energetique ,  Office  National 
d' Etudes  et  de  Recherches  Aerospatiales, 
Chatillon-sous-Bagneux,  France;  test  cases, 
review  of  results 

NASA-Lewis  Research  Center,  Cleveland, 

Ohio;  computation,  test  cases,  review 
of  results 


Air  Force  Aero  Propulsion  Laboratory, 
Wright-Pat terson  Air  Force  Base,  Ohio; 
discussion  and  review  of  designs  for 


Dr.  A.  J.  Wennerstrom 


Dr.  C.  Herbert  Law 

base  and  modified  stages 

General  Electric  Company,  Advanced  Turbo-  Dr.  D.  C.  Wisler 

machinery  Aerodynamics,  Cincinnati,  Ohio;  Dr.  I.  H.  Smith,  Jr. 

review  of  design 

Pratt  6,  Whitney,  East  Hartford,  Connecticut;  H.  We  ingold 
discussion  of  blade  fabrication  problems  A.  W.  Stubner 
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SECTION  VI 


DISCOVERIES,  INVENTIONS  AND  SCIENTIFIC  APPLICATIONS 
No  fundamentally  new  concepts  or  devices  were  developed.  Howevi  , 
Task  IV  involves  some  blade  design  concepts  which  originated  in  AFAPL  and 
may,  after  experimental  evaluation  and  development,  lead  to  improved  multi¬ 
stage  compressor  performance. 
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SECTION  VII 


CONCLUDING  REMARKS 

During  a  two-year  contract  period,  the  program  objectives  supported 
by  USAF/AFOSR  F49620-79-C-0002  have  been  realized.  A  test  installation 
for  generation  of  detailed  flow  field  data  in  curved-mean-line,  rectangular 
cross-section  passages  has  been  built,  instrumented  and  calibrated.  A 
coordinated  set  of  computational  and  experimental  test  cases  is  being 
accumulated  (Task  1) .  New  data  on  unsteady  blade-row  interactions  in  a 
multistage  compressor  environment  have  been  measured,  reduced  and  evaluated. 
These  data  are  expected  to  lead  to  the  introduction  of  mechanisms  for 
unsteady  flow  accounting  in  compressor  design  (Task  II).  A  broad  range 
of  computational  test  cases,  involving  direct  prediction  of  airfoil  cascade 
turning,  loss  and  internal  flow  field  for  compressor  cascade  arrangements, 
was  completed.  These  test  case  results,  when  compared  with  corresponding 
experimental  data  from  linear  cascades,  show  that  an  inviscid-viscous 
interaction  method  is  capable  of  producing  realistic  values.  Predicted 
loss  and  turning  can  be  used  to  study  trends  produced  by  changes  in  most 
of  those  cascade  aerodynamic  variables  that  are  significant  in  current 
compressor  design  (Task  III).  In  a  continuing  project  concerned  with 
evaluation  of  means  for  control  of  secondary  flow  losses  in  multistage 
axial-flow  compressors,  two  two-stage  configurations  were  designed,  a  base 
or  reference  geometry  and  a  secondary  flow  control  stator  geometry.  The 
first  two-stage  blading  set  was  built  and  is  being  installed  in  prepara¬ 
tion  for  the  test  program  (Task  IV).  A  study  was  made  to  recommend  classes 
of  experimental  projects  which  could  contribute  data  useful  to  designers 
and  developers  of  advanced  compressors  and  turbines,  and  which  could  be 
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reasonably  matched  to  the  facility  and  personnel  characteristics  of  uni¬ 
versity  laboratories.  Three  classes  of  experiments  are  suggested  in  this 
report  (Task  V) . 

The  program  in  AERODYNAMICS  OF  ADVANCED  AXIAL- FLOW  TURBOMAOI I NEKY 
has  been  technically  successful.  In  addition  it  has  directed  the  attention 
of  a  significant  number  of  undergraduate  and  graduate  students  toward  the 
objectives  of  current  research  in  a  field  of  aerodynamics  which  lias  not 
received  adequate  emphasis  in  government-university  programs.  Direct 
association  and  cooperation  between  the  aircraft  propulsion  industry, 

USAF  and  NASA  laboratories,  faculty  and  students  will  assisL  in  the  main¬ 
tenance  of  United  States  technological  quality  in  air-breathing  aeronautical 
propulsion. 
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SYMBOLS  AND  NOTATION 


AVDR  ratio  of  blade  row  exit  lo  entrance  axial  velocity  >  density 

product 

b  blade-to-blade  streamsheet  thickness 

blade-to-blade  streamsheet  thickness  at  cas  ade  entrance 

C  pressure  coefficient 

P 

c  blade  chord  length 

D  diffusion  parameter  [Ref.  22] 

HT  total  head 

i  incidence  angle 

M  Mach  number 

M.,  M,  cascade  entrance  Mach  number 
1  1 

m  meridional  distance  measured  from  row  leading  edge 

m  meridional  distance  to  trailing  edge 

PHH  percent  passage  height  from  hub 

Re  Reynolds  number  based  on  chord  length 

Sg  circumferential  distance  between  two  adjacent  stator  blades 

SD  circumferential  distance  between  two  adjacent  rotor  blades 

X  distance  along  chord  line 

Y  circumferential  distance  from  the  reference  meridional  plane 

to  the  measurement  point 

YO  ^  circumferential  distance  from  the  reference  meridional  plane 

to  the  reference  IGV  blade  stacking  axis 

YO^  circumferential  distance  from  the  reference  meridional  plane 

to  the  reference  first  stator  row  blade  stacking  axis 

YO^  circumferential  distance  from  the  reference  meriodional  plane 

to  the  reference  first  rotor  row  blade  stacking  axis 


angle1  of  attack 


ot 

h  local  flow  angle  measured  t rom  axial  direction 

6^  inlet  flow  angle  measured  from  axial  direction 

Y  blade  setting  angle 

k 

6  boundary  layer  displacement  thickness 

k 

9  boundary  layer  momentum  thickness 

0  solidity 

4>  blade  camber 
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APPENDIX  A 


EXPERIMENTAL  DATA  TABULATION  FOR  PASSAGE  TEST  RIG,  TASK  I 

Sample  data  tables,  representative  of  how  data  are  collected  by  the 
data  acquisition  system  and  stored  on  disk,  are  shown  and  discussed.  Three 
data  tables  are  presented: 

(1)  Five-hole  probe  as-measured  traverse  pressures  and  recorded 
temperatures  for  a  given  traverse  plane  and  slide  bar  position. 

(2)  First-pass  reduced  data  based  on  probe  calibration  curves  (Ref.  4). 

(3)  Sidewall  and  endwall  static  pressures  for  given  traverse  plane, 
and  reference  pitot  pressures  at  midpoint  of  entering  boundary 
layer  to  test  section. 

The  first  two  tables  indicate  five-hole  probe  data  for  a  given  "pass" 
of  the  probe  for  various  Z-positions  between  the  endwalls  of  the  passage. 
Probe  position  in  the  passage  is  recorded  by  the  label  in  the  right-most 
column  of  these  tables.  This  label  also  serves  as  file  identification 
for  the  line  of  data  in  the  table  as  stored. 

In  the  third  table,  static  taps  1-16  are  those  in  the  suction  side- 
wall  for  the  given  traverse  plane;  the  taps  are  spaced  at  1-inch  intervals 
with  tap  No.  1  located  0.5  inch  below  the  upper  endwall.  Likewise,  17-32 
are  in  the  pressure  sidewall,  with  No.  17  located  0.5  inch  below  the 
upper  endwall.  ror  static  pressures  along  the  endwall,  obtained  with  the 
static  tap  fixture  installed  in  the  slide  bar,  the  x-position  of  the  tap 
is  recorded  with  the  associated  pressure  measurement. 

Definitions  of  measured  and  reduced  data  terms  are  listed  below  (see 
also  Fig.  A1  for  passage  dimensions,  coordinate  system  and  slide  bar 


numbers) : 


ATM 


PO 

P1,P4,P5 

PBAR 

PITCH 

POS 


PSO 

PT 

Q 

QO 

TBAR 

TPLEN 

VEL 

YAW 


atmospheric  (reference)  pressure,  in  H2O 
plenum  (reference)  pressure,  in  ^0 
five-hole  probe  impact  and  pitch  pressures 
barometric  pressure,  in  Hg  (abs) 

pitch  angle  (deg  determined  from  probe  calibration) 
position  of  five-hole  probe  (three  numbers,  separated  by 
percent  signs,  are  used  indicating  x,  slide  bar  No.,  and 
2,  respectively.  Values  of  x  when  negative  are  prefixed 
by  the  letter  M) 

pitot  static  pressure,  in  ^0  (gage) 

total  pressure,  from  probe  calibration,  in  ^0  (gage) 

dynamic  pressure,  from  probe  calibration 

pitot  dynamic  pressure,  in  H^O 

atmospheric  temperature,  deg  F 

plenum  temperature,  deg  F  or  deg  R 

velocity  V,  ft/s  (see  reduction  formula  below) 

measured  yaw  angle  (determined  by  probe  null)  from  reference 

yaw  for  test  section;  yaw  angle  corrected  for  pitch  (from 

probe  calibration) 


Data  reductions  from  the  tabulated  data  are  as  follows: 
(1)  Total  pressure  coefficient,  C  ,  is  defined  as 


Cpt  =  (pto"Pt)/qo 


where  p  is  total  pressure,  and  p  ,  a  are  reference  pitot  total  and 
t  to  no 

dynamic  pressures,  respectively.  Based  on  table  entries,  then. 
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TRAVERSE 

PLANE  No. 

Y,  INCHES(cm) 

PASSAGE  WIDTH 
INCHES,  (cm) 

1  (INLET) 

-0.500  (-1.27) 

9.00  (22.86) 

2  (15  deg) 

5.498  (13.97) 

9.00  (22.86) 

3  (30  deg) 

10.996  (27.94) 

8.80  (22.35) 

4  (45  deg) 

16.493  (41.89) 

8.48  (21.54) 

5  (60  deg) 

21.991  (55.86) 

8.00  (20.32) 

6  (75  deg) 

27.489  (69.82) 

7.80  (19.81) 

7  (90  deg) 

32.987  (83.79) 

7.80  (19.81) 

8  (DISCHARGE) 

38.487  (87.60) 

7.75  (19.69) 

9  (DISCHARGE) 

43.987  (110.46) 

7.75  (19.69) 

10( DISCHARGE) 

49.487  (125.70) 

7.75  (19.69) 

11 (DISCHARGE) 

54.487  (138.40) 

7.75  (19.69) 

12(01 SCHARGE ) 

60.487  (153.64) 

7.75  (19.69) 

e  A1 .  Schematic  of  test  section  showing  dimensions 
traverse  plane  identification  and  coordinate 
system  convention. 
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T 


(2) 


Static 


C 


pressure 


=  1  -  (PT/QO)  +  (PS/QJ) 
coefficient,  C  ,  is  defined 


as 


C 

ps 


(p  -p  )/q 
s  so  o 


where  p  is  static  pressure  of  the  flow,  and  p  ,  q  are  the  pitot  static 
s  rso  o 

and  dynamic  pressures.  Therefore, 


C  =  (PT/OO)  -  (Q/QO)  -  (PSO/OO) 
ps 


(3)  Flow  velocity  is  determined  from  the  measured  q  and  density  p 
according  to 


With  the  equation  of  state  used  to  express  pg,  the  above  expression 
becomes 


v- 

where  pg,  Tg  are  the  pressure  and  temperature  of  the  flow,  respectively, 
and  R  the  gas  content.  The  pressure  pg  is  in  turn  obtained  from 

Ps  =  [(PT/QO)  -  (Q/QO)](QO)  +  1 3 . 6 (PBAR) ,  in  H?0  (abs) 


Therefore,  the  expression  for  V,  assuming  Tg  =  TPLEN,  and  taking  R  = 
1716  ft^/s^.°R,  becomes 


V  =  58.58 


(Q/QO) (TPLEN) 


[(PT/QO)  -  (Q/QO)  +  13.6 (PBAR/ QO) 


f  t/s 


84 


(4)  Since  the  acoustic  velocity,  c,  is 


c  =  VkRTg  =  tfkR(TPLEN) 
Mach  number,  M  =  V/c,  is  computed  as 


M 


1.20 


(Q/QO) 


f(PT/Q0)  -  (Q/QO)  +  13 . 6 (PBAR/QO) 


in  which  V  from  (3)  above  has  been  used,  and  the  specific  heat  ratio,  k, 
taken  as  1.4. 

(5)  Flow  velocity  components  are  computed  according  to  V  from  (3) 
above  and  the  recorded  pitch  and  corrected  yaw  angles  as 


Vy  =  V/  V 1  +  tan2 (YAW)  +  tan^ (PITCH) 

V  =  V  tan (YAW) 
x  y 

V  =  -  V  tan (PITCH) 

z  y 
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